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ABSTRACT 

We develop models for the self-similar collapse of magnetized isothermal cylin- 
ders. We find solutions for the case of a fluid with a constant toroidal fiux-to-mass 
ratio (r^ = constant) and the case of a fluid with a constant gas to magnetic 
pressure ratio {f3 = constant). In both cases, we find that a low magnetiza- 
tion results in density profiles that behave as p oc at large radii, and at 
high magnetization we find density profiles that behave as p oc r~^. This den- 
sity behaviour is the same as for hydrostatic filamentary structures, suggesting 
that density measurements alone cannot distinguish between hydrostatic and col- 
lapsing filaments — velocity measurements are required. Our solutions show that 
the self-similar radial velocity behaves as oc r during the collapse phase, and 
that unlike collapsing self-similar spheres, there is no subsequent accretion (i.e. 
expansion-wave) phase. We also examine the fragmentation properties of these 
cylinders, and find that in both cases, the presence of a toroidal field acts to 
strengthen the cylinder against fragmentation. Finally, the collapse time scales 
in our models are shorter than the fragmentation time scales. Thus, we antici- 
pate that highly collapsed filaments can form before they are broken into pieces 
by gravitational fragmentation. 

Subject headings: ISM: clouds — ISM: kinematics and dynamics — ISM: mag- 
netic fields — MHD — stars: formation 
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1. INTRODUCTION 

Filamentary structures are commonly seen in giant molecular clouds as sites in which 
are embedded star-forming clumps. These structures are threaded by magnetic fields (Heiles 
et al. 1993; Matthews & Wilson 2000; Falgarone et al. 2001), which can support the filaments 
against gravitational collapse. While the mechanism by which filaments form is not well 
understood, it has been suggested that these filamentary structures occur as the result of the 
fragmentation of the parent molecular cloud through turbulent motions (Klessen & Burkert 
2000; Ostriker et al. 2001). Recent theoretical work (Fiege & Piidritz 2000a,b) suggests that 
helical magnetic fields may play an important role in determining the structure of a filament. 

The study of the collapse of a spherical pre-stellar core has been greatly benefitted by 
the study of one-dimensional models (Larson 1969; Penston 1969; Shu 1977; Foster & Cheva- 
lier 1993), which have become very useful in identifying the physical origin of a particular 
structural phenomenon. By leaving out most of the complicated and analytically intractable 
physics which may occur in the full three-dimensional system, these models allow us to iden- 
tify the physical processes which lead to the observed structure, and can be compared to 
both numerical simulations and observational data of molecular cloud cores. The spherical 
self-similar models feature an initial gravitational collapse phase that ends with the appear- 
ance of a singularity in the density field. This is followed by an accretion phase, wherein 
an expansion wave propagates outward from the newly formed core enabling the infall of 
the remaining envelope. Similar analyses of the process of filament formation and evolution 
have not been as well studied, despite the frequency of filamentary structures forming in 
numerical simulations (Porter et al. 1994; Klessen & Burkert 2000; Ostriker et al. 2001), and 
the abundance of filaments in observations of molecular clouds (e.g. Falgarone et al. 2001). 

Stodolkiewicz (1963) found solutions for the equilibrium structure of an isothermal cylin- 
der, for a fluid where the ratio of gas pressure to magnetic pressure is constant. Changing 
the value of this ratio has the effect of changing the scale radius of the solution, which falls 
off as (r/ro)~^. Ostriker (1964) found identical solutions for the unmagnetized isothermal 
cylinder. 

Miyama et al. (1987) derived a set of self-similar solutions for an unmagnetized isother- 
mal cylinder. These solutions have an infall velocity which is proportional to the distance 
from the axis, and a density structure that has the same form as the unmagnetized equi- 
librium fllaments of Stodolkiewicz (1963) and Ostriker (1964), but with a scale radius that 
decreases with time and a central density that increases with time. 

Recently, Hennebelle (2002) explored a set of self-similar solutions for a rotating, mag- 
netized fllament which may undergo collapse in the axial direction, in addition to radial 
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coUapse. This set of solutions include the equations which we will derive in the next section, 
although Hennebelle (2002) does not explore the case where there are only purely radial 
motions. This is the focus of our analysis. 

The goal of this paper is to extend the calculations of the self-similar collapse of an 
isothermal cylinder to include magnetic and rotational effects. We show that one has to 

choose a relationship between the magnetic field strength and the density. When this is done, 
we find that the density and magnetic field structures of collapsing cylinders is identical to 
that of the corresponding hydrostatic structure, with only the scale radius and normalization 
of the quantities changing in time. The important consequence of this for observations is that 
only careful large-scale velocity measurements can discriminate between static and evolving 
filamentary structure, unlike the situation for accreting spheres. We derive and solve our set 
of self-similar equations in Section 2. These solutions will be discussed in Section 3. 



2. THE SELF-SIMILAR COLLAPSE EQUATIONS 



We begin our analysis with the equations of ideal magnetohydrodynamics for a self- 
gravitating gas: 



dp 
dt 



+ V-(pv) = 

- = Vx(vxB) 



p— p(v • V) V 



dt 



-VP + 



(B ■ V) B - - V (B^) 
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(1) 
(2) 

(3) 
(4) 



In these equations, p is the density, v is the velocity, B is the magnetic field, P is the thermal 
pressure, for which we adopt the isothermal equation of state P — c^p (with a constant speed 
of sound Cg), and $ is the gravitational potential. 

The geometry of our filaments leads us to use a cylindrical coordinate system to solve 
the equations. Our cylindrical symmetry causes all d/d(f) and d/dz derivatives to be zero, so 
that the resulting fields are functions of (r, t) only. This geometry allows us to replace the 
gravitational potential by a mass per unit length m, defined by 



dm ^ 
— = 27rrp 
or 



(5) 



We can make the following substitutions, to change our equations from functions of (r, t) 
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to separable functions of {s,t), for the self-similar variable s — r/{cst). These solutions will 
represent collapsing motions for s < (and hence t < 0), and expanding motions for s > 0. 



P{r,t) = (6) 
m{r-,i) = —^m{s) (7) 

Gr 

v(r,f) = c.v(s) (8) 

^('■•') = y^T«W 

where p{s) etc. arc dimcnsionlcss functions of the self-similar variable s alone. Substituting 
these into Equations (l)-(5) (see Appendix A for details), 

Vr = S (10) 

dm , , 

d 



= (12) 
= (13) 
4§ = (14) 



s ds s s ds \ / ds 

Br d f ~ \ 



Before we proceed to solve Equations (10)-(17), we will first recap the solutions to 
these equations in the case of B = 0. Equation (11) remains unchanged, Equation (17) 
shows that Vip = 0, and Equation (16) reduces to d(ln p)/d(lns) -|- 4m = 0. The solution to 
this, combined with Equation (11), is p = 2C (1 + C^^) and m = (1 + C^^) , where C 
is a constant of integration. This solution was first found by Miyama et al. (1987). It has 
two features worth highlighting: 



1. It has the same radial density structure as the equilibrium solution of Stodolkiewicz 
(1963) and Ostriker (1964). 
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2. The total enclosed mass per unit length (i.e. as |s| — > oo) is the critical mass per unit 
length. 

The asymptotic density behaviour of this solution tends to at large values of s; we 
will show later that magnetic solutions tend to have this behaviour only in the limit of low 
magnetization. 

There are two classes of solutions to Equations (10)-(17) when the effects of magnetic 
fields are included, which we can separate through the use of Equation (14): 

1. Br = 0: In this case, we can show from Equation (15) that Bz = and from Equation 
(17) that = as well. This solution also has the consequence that our system of 
equations is underspecified, as B^ only appears in the radial momentum equation. We 
will discuss the solutions to this below. 

2. dBz/ds = 0: If our Bz field is constant, then from Equation (12) wc can see that we 
have a radial field Br = Bro/s which is singular along the cylinder axis. 

The latter solution is difficult to justify physically; even though the field is still divergence- 
free, the field lines are not closed because of the linear singularity at the axis of symmetry. 
Henceforth we discard this possibility and focus our attention on Case 1. 

We now proceed to solve the remaining set of equations. With the condition that Br = 0, 
we are left with 

(18) 

sp (19) 

f + ^ + ^±(sB,) (20) 
as s s as \ / 

We can integrate Equation (20) to obtain the following relationship between p, m and 

B^: 

a = (^2p + S|) + 4m (m - 1) (21) 

We can constrain the value of the constant a by examining the behaviour of the system on 
the axis. Since limm = 0, it follows that in the limit s —>■ 0, lims^(2p + B^) = a. 

The first term in this result can be evaluated by using Equation (19). We can express 
the density asp = A|s|"ass^O, as for sufficiently small s, one term in the expansion of p 



Vr = 

dm 
ds 

= 
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will dominate the other terms. Prom Equation (19), we can write the mass per unit length 
for this density behaviour as m = A jj*' |s|'^"'"^d|s|. For n < —2, fh will not be finite in the 
limit s — > 0, so with the constraint n > 2, the first term in Equation (21) has the limit 
lim ps^ — 0. 

s->-0 

We can constrain the magnetic field behaviour in a similar manner, by examining the 
toroidal magnetic flux, ^. We use the formulation of the magnetic flux of Fiege & Pudritz 
(2000a), such that ^ = J iJ^da, where da = dr dz is the area element through which the 
magnetic field passes (sec Figure 1). If we consider the fiux per unit length ip = d^/dz, and 
write the magnetic field as B^j, = C|s|"' as s 0, we find that if) = C |s|"'ds. For n' < — 1, 
ijj is not finite, so we restrict n' > —1. We thus establish that a = in Equation (21). 

The properties of the system can be examined if we solve for the mass per unit length 
in Equation (21). 



1 1 



= ^±-Jl-s^{2~p + Bl) (22) 



We can identify two branches to this solution. The minus branch m_ < 1/2 corresponds 
to the minus sign in Equation (22); the plus solution m+ > 1/2 is the other branch. The 
two branches meet at the half-mass point given by s^^ 2p(si/2) + -B|(si/2) = 1- Since 

m(0) = 0, we know that the critical quantity gcrit = s'^i'^P^ B"^) — >^ as s ^ 0. Thus, we can 
determine that for every solution, the quantity gcrit must increase from — > 1, then begin 
to decrease again. Thus, the minus branch begins at 7ft_(0) = 0, then increases to 1/2, then 
decreases again; the plus branch starts at m+(0) = 1, decreases to 1/2, then increase again. 

Since the plus branch is always greater than 1/2, it cannot represent the solution at 
s = 0; similarly, at values of \s\ > \si/2\, the mass per unit length must continue to increase 
(as a decreasing mass per unit length implies a negative density according to Equation (19), 
which is unphysical), and so the minus solution is unphysical in this regime. Thus the 
physical solution must begin on the minus branch at low s, then at the half-mass point 
switch to the plus branch. The maximum value of fh must occur when p — 0. If this occurs 
at a finite radius Smax, then the total mass per unit length will be 



^max = ^±^V^1-Smax^|(w) (23) 

where the plus or minus sign is determined by whether jsmaxl > \si/2\ or Ismaxl < \si/2\, 
respectively. This rfimax is identical to the reduced critical mass per unit length mmag found 
by Fiege & Pudritz (2000a) from the virial theorem, for a filament with a toroidal magnetic 
field. 
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Fig. 1 



, — Area element for the calculation of the toroidal magnetic flux. 
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Because of the indeterminacy in this solution as noted above, we need to specify a form 
for in order to make further progress. We will perform this calculation for two plausible 
magnetic field configurations which commonly appear in star formation literature: one where 
the ratio of the thermal to magnetic pressure f3 = 2fic'j.p/B'^ = 2p/B'^ is constant, and one 
where the toroidal flux-to-mass ratio = B^/{\s\p) is constant (Fiege & Pudritz 2000a). 
We will discuss these two cases in turn. 



3. RESULTS 

There are two models for the magnetization of the gas in molecular clouds which are 
used extensively in the literature: one in which the magnetic pressure is proportional to the 
gas pressure, and one in which the magnetic flux through a fluid clement is proportional to 
the mass of that element. As we have the freedom in our set of equations to specify one of 
our variables, we choose to constrain the magnetic field through one of these relations. 

As we show in the Appendix, the structure of the self-similar collapse equations (19)- (20) 
is identical in form to the equations for the hydrostatic solutions; the density and magnetic 
fields have identical radial structures. The only difference between the static and collapse 
solutions is that the collapse solutions are time-dependent, and have an infall velocity. It 
is remarkable that the radial velocity terms in the momentum, induction and continuity 
equations exactly cancels out the time-dependent terms. A very different behaviour is seen 
in spherical symmetry, for which the collapse solutions have a different density structure than 
the static solutions. This underlies the importance that the geometry has in determining 
the spatial structure of the collapsing filament. 

In self-similar cylindrical coordinates, the continuity equation (1) restricts the radial 
velocity to a single solution, Vr — r/t. Thus, this analysis can only represent a collapsing 
solution for t < 0. Furthermore, this implies that there is no equivalent to the accretion 
solution found in the self-similar collapse of a spherical isothermal fluid by Larson (1969), 
Penston (1969), and Shu (1977). The collapse of filaments, therefore, is very different than 
the collapse of spheres. 

3.1. Case lA: Constant P 

A common description for magnetic fields in molecular clouds is a constant ratio of gas 
pressure to magnetic pressure, j3, such that the magnetic field is B^ = {2j3~^p)^^^. We can 
use this in Equation (22) to write 
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Fig. 2. — Structure of the (3 — constant gas. Panel (a) shows the density structure; (b) shows 
the behaviour of the mass per unit length, m; (c) displays a log-log plot of the density; and 
(d) displays the toroidal magnetic field strength B^. The thick solid line in (c) represents the 
curve p oc s^^; the thick dashed line represents the curve p oc s^"^ (the asymptotic solution 
to the non-magnetic case). 
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rn = l±^^l-2s^{l + p-^)p (24) 

We can differentiate this and use Equation (19) to solve for p(s), with the result 

sech^ [(ln|s|)/(l + /3-^)] 
'^'^ ^ 2sHl + p-^) ^''^ 

rh{s) = ^{l + tanh[(ln|s|)/(l + /?-^)]} (26) 

~ _ sech[(ln|.|)/(l + r^)] 

""'^'^ - — i^wm — ^ ^ 

This solution displays a number of the properties mentioned previously. The half- mass point 
occurs when |s| = 1; the value of 5^^2(2^(51/2) + B'^isi/'i)) — sech^(O) = 1, as expected. This 
is also the sonic point in the fluid (that is, where \vr \ — 1). The limit lim (sB^) — 0, so the 

|s|— >oo 

total mass per unit length ought to be equal to the critical value rrimax = 1; from Equation 
(26) we can see that this is indeed the case. The solutions to Equations (25)-(27) are plotted 
in Figure 2. The behaviour of the quantity g'crit — {^P + is plotted in Figure 3a, and 
displays exactly the behaviour discussed in Section (2). The branches of fh are shown in 
Figure 3b. 

Since most measurements of molecular clouds have difficulty resolving the interior re- 
gions of cores and filaments, we examine the behaviour of these solutions at large s (equivalent 
to large r). We can do this by expanding the sech^ term into exponentials. The net result is 

r, 2/3 + 1 

poc |s| ^ (28) 

for |s| ^ 1. Equation (28) shows that for /3 ^ 1 (corresponding to weak magnetic fields), 
p oc s~^. For /3 <^ 1 (strong magnetic fields), p oc Note that in the limit s — > 00 

(equivalent to t — > 0), p ^ for all r 7^ 0. This simply shows that there cannot be an 
accretion solution for t > — there is no material left to be accreted. 

We can similarly write an expression for the density behaviour as s — > 0: 

p oc (29) 

For a finite value of (3, the density will be singular at the axis. This is in contrast to the 
result of Stodolkiewicz (1963), who used an arbitrary boundary condition of p{r = 0) = 
constant. We beheve our result to be physically more self-consistent, as it satisfies the 
physical constraints of zero mass per unit length and toroidal flux at the origin (the solution 
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of Stodolkiewicz (1963) cannot simultaneously satisfy the hydrodynamic equations and the 
restriction of m(0) = 0, for a purely toroidal field). Our solution shows that the density in a 
(3 = constant model will be singular on the axis. A singularity in the spherical solutions can 
be rationalized as an approximation to the effects of a protostar in the center of the cloud. 
In the cylindrical case, even though collapse as a singular filament occurs, this structure 
itself will undergo fragmentation in a finite amount of time, as is shown in Section 3.3. Our 
treatment of the P = constant case shows that it may be of hmited astrophysical interest. 

The toroidal field contributes both a magnetic pressure force and a magnetic tension 
force to the momentum equations. It is instructive to compare the magnitudes of these two 
forces. The magnetic pressure is (1/2) d(i?|)/ds, and the magnetic tension force is —B'^/s. 
The ratio is thus 

R = 1 _l dHBl) .30^ 

2(B-V)B 2 dln\s\ ^ ' 

For our (5 = constant fluid, R = (1/2) d(lnp)/d(ln |s|) = 2{2m + (3-^). The magnetic 
pressure force will dominate for \R\ > 1, which corresponds to m > 1/4 — (l/2)/3~^; note 
that this will always be the case for a fluid with /? < 2, as then the right-hand side will 
be negative. Strongly magnetized filaments are therefore supported by the pressure forces 
associated with their toroidal fields. 



3.2. Case IB: Constant Flux-to-Mass Ratio 

Another form for the magnetic field configuration of a filament that has been proposed 
(Fiege & Pudritz 2000a) is one wherein the magnetic fiux through a volume is proportional 
to the mass in the volume, such that the fiux-to-mass ratio is constant everywhere {B^ = 
F(^|s|p). With this assumption, we can numerically integrate Equation (20) using a Runge- 
Kutta integrator to obtain a set of solutions. We plot these solutions for several values of F^ 
in Figure 4. As for the solutions for the /? — constant plasma, the solutions for F<^ = constant 
exhibit the properties discussed in Section (2). However, we also note some differences in the 
behaviour of these two different systems. We plot the critical quantity gcrit for F<^ = constant 
in Figure 3c, which can be directly compared to the g'crit of the (3 — constant gas, plotted 
in the panel above it. While the half-mass point for the f3 = constant case did not depend 
on the value of the parameter f3, for the F^ = constant solution the half-mass point does 
depend on the value of F^. The density and magnetic field for the F<^ = constant solution 
are also non-singular, as can be seen in Figure 4. 

At low values of magnetization (i.e. low F^), the density distribution is very similar 




Fig. 3. — Behaviour of the critical quantity g'crit and the mass per unit length rh. The upper 
panels show the behaviour of the [3 = constant gas; the lower panels are for the = constant 
case. The panels on the left display the behaviour of the critical quantity s^(2p + i?|); those 
on the right show the upper branches (dotted lines) and lower branches (solid lines) of the 
mass per unit length. The heavy lines mark the allowable branches. The different lines in 
(a) describe the same models as the lines in Figure 2; the different lines in (c) describe the 
same models as the lines in Figure 4. 
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to that predicted by Miyama et al. (1987) for the self-similar collapse of an unmagnetized 
filament (which in turn has a the same shape as the unmagnetized hydrostatic filament). As 
the magnetic flux increases, the density distribution at large |s| tends towards p oc rather 
than the steeper s~'^ for weak magnetic fields. This can be shown explicitly by isolating dp/ ds 
in Equation (20): 

dp ^ _f2p\2m + ryp 

ds \s J 1 + ryp ^ ' 

For -C 1, we see that dln{p)/dln\s\ = —4m; since m ^ 1 as |s| — >• oo, p oc in 
the case of low magnetization. For » 1, we see that dln{p)/dlB. \s\ = —2, implying that 
p oc as discussed above. These results and those in Section 3.1 show that the asymptotic 
behaviour at large radii for very weak or very strong magnetic fields is independent of the 
details of how these filaments are magnetized. 

As for the constant-/? solution, as s — > oo, p — > 0. Thus, for solutions with a constant 
F^, there cannot be an accretion solution for t > 0, either. 



3.3. Fragmentation 

Stodolkiewicz (1963) and Fiege & Pudritz (2000b) derived criteria for the fragmenta- 
tion of a magnetized isothermal cylinder, for each of the forms of the magnetization that we 
have used. Since our collapse solution moves towards a line singularity at t=0, which is un- 
physical, there exists the possibility that the cylinder will become unstable to fragmentation 
instabilities as the density and magnetic field strengths change with time. 

Fiege & Pudritz (2000b) found that for a filament with constant toroidal fiux-to-mass 
ratio, the fragmentation timescale is 

where w^^^, is a frequency of fragmentation which depends on the value of F<^. Fiege & 
Pudritz (2000b) found that -cul^^x < 0.0581. 

As the radial density structure is the same in both the self-similar case and the hydro- 
static case, we can compare the fragmentation timescale to our collapse timescale, to see 
which process should dominate. The collapse timescale is 

1 / n \ -1/2 
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Fig. 4. — Structure of the — constant gas. Panel (a) shows the density structure; (b) 
shows the behaviour of the mass per unit length m; (c) displays a log-log plot of the density; 
and (d) displays the toroidal magnetic field strength i?^. The thick solid line in (c) represents 
the curve p oc s~^; the thick dashed line represents the curve p oc s^"^ (the asymptotic solution 
to the non-magnetic case). 
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and so the ratio of the fragmentation to collapse timescales is 

, ,2 \ -1/2 



^/-p^4_75j_:ii^| (34) 



Tcoi V 0-01 

This ratio has a minimum at F,^ = of Tfrag/Tcoi = 1-97, which increases as increases. 
Thus, the toroidal magnetic field acts to prevent the cloud from fragmenting. In the case 
of low magnetization, the collapse and fragmentation timescales are approximately equal, 
suggesting that the cloud may begin to collapse along the filament axis, in addition to its 
radial collapse. 

We can also find the ratio of the critical wavelength for fragmentation to the radial 
collapse scale. Fiege & Pudritz (2000b) found that the fragmentation wavelength for a 
constant toroidal flux-to-mass ratio is 

where kmax is a dimensionless wavenumber corresponding to the maximum instability, and 
has a maximum value of 0.462 at F^ = 0. An estimate of the radial length scale can be set 
by setting |s| = 1: 

^rad — CgSt 



c 



From this, we get a ratio of fragmentation to radial wavelengths of 

^ ^ 14.43 ' (37) 

Kad \ 0.2 / 

At = 0, this has its minimum value of 6.25 (3.41 if we use the critical wavelength rather 
than the maximum wavelength). Again, this suggests that toroidal fields support the cylinder 
against collapse. 

We can do a similar analysis for the (3 = constant case. The fragmentation wavelength 
found by Stodolkiewicz (1963) for a = constant fluid is 

while the radial wavelength is still given by Equation (36), for a ratio of 

^ - 5.563./l^ (39) 

Kad V P 
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Again, this has a minimum value when the magnetization is zero (/3 = oo); the presence of 
a toroidal field serves to stabilize the cylinder against fragmentation. 

In addition to the axisymmetric instabilities discussed above, there exists the possibility 
that non-axisymmetric instabilities (e.g. kink modes) could develop. The presence of such 
instabilities would destroy the axisymmetric nature of our models, a problem that cannot 
be addressed with a self-similar formalism. We can, however, estimate the likelihood of such 
modes by comparing their growth rates with that of the axisymmetric modes. Nakamura 
et al. (1993) has investigated the kink mode instabilities of a hydrostatic molecular cloud 
with helical magnetic fields described by a constant i3. They found that the axisymmetric in- 
stabifitics, such as those described above, have faster growth rates than all non-axisymmetric 
modes. Thus, axisymmetric modes will dominate the fragmentation, and non-axisymmetric 
modes are not an important effect in the evolution of these filaments on the collapse time 
scales that we are considering. 



4. DISCUSSION AND CONCLUSIONS 

We have determined the structure of a collapsing self-similar magnetized isothermal 
cylinder for the two most common magnetic field configurations found in the literature. 
One of the major conclusions of our study is that the cylindrical radial structure for the 
gas density is indistinguishable from the hydrostatic structure with the same magnetic field 
configuration. Consequently, one would need to examine the profiles of molecular lines for 
infall signatures in order to differentiate collapsing and equilibrium structures. Because the 
high-density regions that emit most of the thermal radiation will have low velocities, it may 
be more difficult to detect a characteristic line profile in comparison with the situation for 
accreting spherical cores that can be identified by an inverse P-Cygni profile. 

The self-similar density profile at large radii from the axis follows a power law. For a 
cylinder with a weak magnetic field, we find that at any particular time during the collapse, 
p oc r^"^. By way of contrast, for a cylinder with a strong magnetic field, we find that 
p oc r~^. These results are independent of the model for the magnetic field strength that we 
use {P = constant or = constant). 

Our solutions require that there be no motions along the axis, and that the magnetic field 
is purely toroidal. In a fully dynamic, three-dimensional collapse, we would not necessarily 
expect this to be the case. However, if the axial motions and poloidal fields are sufficiently 
weak, we could reasonably expect our solutions to be a good approximation. Furthermore, 
the static models of Fiege & Pudritz (2000a) indicate that for models consisting of both 
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poloidal and toroidal fields, the toroidal magnetic field strongly dominates in the outer 
regions of a filament, suggesting that our solutions may be applicable to the envelopes of 
filaments threaded with helical fields. 

The velocity behaviour of our solutions bears a close resemblence to the solution found 
by Penston (1969) for the self-similar collapse of an isothermal sphere, in that the infall 
velocity Vr —>■ as r ^ 0. However, due to the simple nature of the continuity equation, we 
only find the collapse solution; there is no analogue to the Shu (1977) accretion-expansion 
wave solution that follows the formation of a singular core. It is also the continuity equation 
which requires f oc r at all radii, rather than having an asymptotic infall velocity as found 
by Penston (1969). Clearly, arbitrarily high velocities at large radii are not physical; there 
needs to be a mechanism to truncate the models at finite r. Pressure truncation, which we 
have not included in this analysis, could satisfy this condition. 

We find, in agreement with Stodolkiewicz (1963) and Fiege & Pudritz (2000b), that 
the presence of a toroidal field stabilizes the cylinder against fragmentation. Even in the 
absence of a field, however, collapse motions should dominate over fragmentation. We thus 
have shown that self-similar filaments should remain relatively unfragmented as they collapse. 
This mechanism would allow filamentary structures to survive over long periods of time while 
retaining their filamentary character. If this were not true, we would not expect to see many 
filaments because they would fragment into cores too quickly. 

Although large scale velocity fields have been observed to be associated with individual 
cloud cores (e.g. Mardones et al. 1997; Lee et al. 1999), we emphasize that the flows 
predicted here arc on a yet larger scale associated with filament formation. Simulations of 
hydromagnetic turbulence show that filaments can form from shocks (Ostriker et al. 1999, 
2001). We plan to investigate the possible link between our analytic models and time- 
dependent simulations of magnetized clouds. 

We thank an anonymous referee for useful comments on the manuscript. D.A.T. is 
supported by an Ontario Graduate Scholarship. R.E.P. is supported by the Natural Sciences 
and Engineering Research Council of Canada. 

A. THE STRUCTURE OF THE SELF-SIMILAR AND HYDROSTATIC 

EQUATIONS 

We can show explicitly that the structure of a filament that is collapsing self-similarly 
is identical to that of a hydrostatic filament with the same constitutive equations for the 
magnetic field. In our formalism, we replace the gravitational potential with the mass per 
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unit length (Equation 5) 
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2Gm 



(Al) 



We write out Equations (l)-(4) in component form, and use the fact that we are considering 
only the axisymmetric and axially symmetric case (i.e. functions depend only on the radius 
r, and the time t). 
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If we make the substitutions given in Equations (6)- (9), we obtain the self-similar struc- 
ture equations (10)-(17). If we instead consider the equilibrium solutions (so that v = 0, 
d/dt = 0), we get 
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From the condition that V ■ B = 0, we require Bj. = 0, for the same reasons we had in 
Section 2 for Equation (12). We now convert to dimensionlcss quantities, using (in analogy 
with Equations 6-9) 

1 



P = 



(A13) 
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(A14) 
(A15) 

(A16) 
(A17) 



This results in 



dm 
dx 



xp (A18) 

^ dp ^ Amp ^ d ^ ^ ^^ _^ dB^ (A19) 

dx X 

Comparing Equations (A18) and (A19) with Equations (19) and (20), we see that in the 
cylindrically symmetric geometry, the self-similar collapse equations and hydrostatic equi- 
librium equations will always give the same structure for the density (with the exception 
that a Bz field in the hydrostatic case is allowed) , given a constitutive relation between the 
toroidal magnetic field and the density. 
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